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In a point contact NS junction, perfect Andreev reflection occurs over a range of voltages equal 
to the superconducting energy gap, producing an excess current of /oxc ~ (4/3)(2eA//i). If the 
superconductor has a finite width, rather than the infinite width of the point contact, one cannot 
neglect superfluid flow inside the superconducting contact. The energy range available for perfect 
Andreev reflections then becomes larger than the superconducting gap, since superfluid flow alters 
the dispersion relation inside the finite width superconductor. We find a maximum excess current 
of approximately (7/3)(2eA//i) when the width of the superconductor is approximately 7/3 times 
the width of the normal metal. 



I. INTRODUCTION 



At a normal metal - superconductor (NS) junction, 
electrons incident from the normal metal can be scat- 
tered into time reversed electrons (holes) by the pairing 
potential. -This conversion process is known as Andreev 
reflection.EJ When the NS junction carries an electrical 
current, Andreev reflection is accompanied by the conver- 
sion of normal current to supercurrent fl'u This supercur- 
rent flow modi^p. the current- voltage relaticin in NS and 
NSN|-|iunctionsCElj, superconducting wiresD'El, SNS junc- 
tionscl, and NS umctions with a supercurrent parallel to 
the NS interfacdlj. In this paper we consider the current- 
voltage relation for NS junctions having a supercurrent 
flow perpendicular to the NS junction. 
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FIG. 1. Different types of NS Junctions, (a) NS point con- 
tact, (b) N-Narrow S- S, (c) N - Wider S- S. The wire width 
determines the number of conducting modes in the narrow 
segments {Mm and Ms). The superfluid flow velocity Vs can- 
not be neglected when Mm — Ms- 

The superfluid flow present for the point contact NS 
junction in Fig. ^a) will have little effect on its I-V 
characteristic, since the current density inside the wide 
superconductor approaches zero. However, for the NS 



junctions shown in Fig. |l|(b)-(c) , the number of con- 
ducting modes in the superconducting wire {Ms) is com- 
parable to to the number of conducting modes in the 
normal wire {Mn). Since the current density is not zero 
inside the superconductor, one cannot neglect the effect 
of a superfluid flow on the I-V characteristics of the NS 
junctions shown in Fig. |l|(b)-(c). Since the superfluid 
flow strongly modifies the dispersion relationship of the 
superconductor when the ratio of the number of conduct- 
ing modes a = Ms/Mn is of order one, including such 
a superfluid flow will influence the I-V characteristic of 
NS junctions. Since the number of conducting modes 
is roughly proportional to the width of the conductor, 
namely a ~ Ws/Wn, we can vary the ratio of conduct- 
ing modes by varying the width of the superconductor 
Ws. 

Blonder, Tinkham and Klapwijk (BTK) showed that 
the point contact NS junction in Fig. |l|(a) carries a lareer 
current than a normal metal point contact junctiontll. 
This 'excess' current of is due to the presence of Andreev 
reflection, and has the value /oxc = (4/3)(2eA//i) in a 
ballistic point-contact NS junction. By varying the width 
of the superconducting wire forming the NS junction in 
Fig. |l|(c) , we determine in this paper how the excess cur- 
rent varies as a function of the transverse mode ratio a. 
Using both an independent band model, and a second 
model which includes scattering between different lateral 
modes at the NS interface, we show the excess current 
in NS junction can be much larger than the BTK result. 
This enhancement of the excess |-current over the BTK 
value has also been noted in Ref.El. 

More Andreev reflections can occur at higher volt- 
ages when a supercurrent flows perpendicular to the NS 
interface, accounting for an excess current larger than 
the BTK result. If the superconductor is narrower than 
about Ws < (7/3)M^Ar, the narrow region of the super- 
conductor exceeds its critical current before allowing the 
maximum number of Andreev reflections. If the super- 
conductor is much wider than Ws >> {7/3)Wn, there is 
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too much geometrical dilution of the supercurrent for a 
significant superfluid flow to develop inside the narrow- 
est region of the superconductor. We find a maximum 
excess current when the width of the superconductor is 
approximately Ws — (7/3)M^Ar. 



II. WHY SUPERFLUID FLOW INCREASES THE 
EXCESS CURRENT 



In this section we give the simplest physical model 
which illustrates why the excess current in NS junctions 
can be larger than the point contact limit. To make our 
physical points, we construct a crude two fluid model, 
which does not obey electrical current conservation at 
every point in space. A fully self-consistent solution 
of the Bogoliubov-deGennes equations automatically en- 
sures current conservation, eliminating the need for this 
ad-hoc two-fluid model, but requires more computational 
effort. The two-fluid model we develop only guarentees 
current conservation at the terminal contacts. A fully 
self-consistent solution of the BdG equations in a sin- 
gle mode junction (a = 1), done in section [II, 



gives 

the same results for the excess current as the two-fluid 
model. Viewing the electrical conduction in terms of this 
two fluid model, therefore, allows us to obtain a value of 
the superfluid velocity Vg inside the superconducting con- 
tacts for each value of the bias voltage across the NS or 
NIS junction. We then use this value of the flow velocity 
Vs to compute Andreev and normal reflection probabil- 
ities for each value of the voltage, and thus obtain the 
electrical current in a globally self-consistent manner. 



A. Two Fluid Model 



Figure g shows the energy band diagram of an NS 
junction when the superconductor carries a finite super- 
currentij. In a transmission formalism, one must com- 
pute the electrical current operator for all incident quasi- 
particle states, and add them to obtain the total cur- 
rent. Ref£3 evaluates the electrical current operator on 
the normal side of the NS junction in terms of parti- 
cle current transmission and reflection probabilities. The 
derivation in section 3 of RefJlj is valid for a multipjle 
moded NS junction subject to a superfluid flw- Ref.tHi 
recovers the well-known BTK current formulaU, namely 
(for zero temperature) 

I = (2e//l)Mjv / [1 ~ TNe.NeiE) + TNh,Ne[E)]dE. 

Jo 

(1) 

In Eq. ([^) we use the notation of Ref.lli, where the Tns.np 
particle current reflection probabilities from the incident 



channel {N jS) to reflected channel {N5). The indices 
/3,S — e or h for electron-like or hole-like quasi-particles. 
We must compute the both the normal Tj^e.Ne — Rn 
and Andreev Tj^h.Ne — Ra reflection probabilities when 
the superconducting contact is subject to a superfluid 
flow. For simplicity we have taken the modes in the 
normal conductor to be both independent and to carry 
identical electrical currents. 
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FIG. 2. Energy band diagram of an NS junction sub- 
ject to a superfluid flow vs > 0. The shifted energy bands 
in the superconductor cause Andreev reflection to occur at 
higher energies than without superfluid flow. The contacts 
inject electron-like (solid dots) and hole-like (open circles) 
quasi-particles as shown. 

In order for the superconductor to carry a finite super- 
current, the form of the order parameter A(2;) inside the 
superconductor must be generalized to A(x) = |A|e^*'^. 
The superfluid velocity is Vs = hq/m. If we consider 
solutions accurate within the Andreev approximation 
(|£'| ^ /i), we can approximate the dispersion relatio 
near the Fermi level as being rigidly shifted in energy a: 
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(2) 



Here ^Qd/'m — \A\/hkp is the Landau depairing ve- 
locity.E3 The electrical current is given by Ic — enVgMs 
with n = 2{kp/TT) the electron density per mode and 
Ms the number of (equivalent, fer simplicity) conduct- 
ing modes in the superconductoiCl. We can rewrite the 
supercurrent Ic in terms the depairing velocity as 



Ic ^ {Ae\A\/h){v,/vd)Ms 



(3) 



At zero temperature, the critical current phase bound- 
ary occurs when Vs = Vd- We must therefore maintain 
Vs < Vd to preserve the superconducting order parameter 
|A|^0. 

Inside the superconductor, the electrical current is of- 
ten argued to be com poso id of a 'quasi-particle' and 'con- 
densate' contributionllilliil. To break Eq. dl^) down into 
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these-jtwo contributions to the current we use the sum 
ruleB 

1 — TNe,Ne{E) + TNh,Ne{E) + Tse.Ne{E) + Tsh,Ne{E). 

(4) 

Equation (Q) states that the normal {Tse.Ne — TV) and 
Andreev {Tsh.Ne = Ta) particle current transmission co- 
efficients into the superconductor conserve the total num- 
ber of quasi-particles. Combining Eqs. (|l|) and (||), the 
electrical current inside the superconductor at zero tem- 
perature is 



QP 



I A. 



(5) 



We identify the portion of the electrical current due to 
'quasi-particlc' injection as 



eV 



i2e/h)MN [Tse.Ne{E)^Tsh.Ne{E)]dE, (6) 



and the 'Andreev' portion of the current Ia as 

Ia = {2e/h)MN / 2[TNh,Ne{E) +Tsh,Ne{E)]dE. (7) 
Jo 

The 'Andreev' current in Eq. (|^) equals twice the sum of 
all Andreev processes. 

Equations give the same current as the BTK 

expression from Eq. (|^), since both are simply the cur- 
rent operator evaluated inside the normal metal. The 
key physical element in our 'two-fluid' model is that we 
require that both Eqs. (||) and (0) for the 'condensate' 
current Ic and the 'Andreev' current I a must be equal. 
We examine this assumption more rigorously in the Ap- 
pendix, by evaluating the electrical current operator in- 
side the superconductor. In this section we make a plau- 
sibility argument for equating Ia from Eq. with the 
'condensate current' Ic from Eq. (||), or the 'current of 
Cooper pairs'. 

Eq. (0) expresses the condensate current in terms of 
probabilities for Andreev reflection and Andreev transmi- 
sison of an electron incident from the normal metal. Con- 
versely, Eq. (H) expresses the condensate current in terms 
of the superfluid velocity. That an incident electron and a 
reflected hole on the normal metal side requires a Cooper 
pair to move off into the superconductor (to-preserve 
electrical charge conservation) is well knowntll. Simi- 
larly, Andreev transmission of an electron also requires 
a Cooper pair flow inside the superconductor. Andreev 
reflection and Andreev transmission of electrons incident 
from the normal metal require Cooper pairs flow away 
from the NS interface (in the same direction) for both 
processes. Equation (Q) embodies this physical reason- 
ing. Similarly, identifying Eq. ^ as the 'quasi-particle' 
contribution to the electrical current is also quite natu- 
ral, being proportional to the electrical current operator 
for an electron injected from the normal metal (evaluated 
on the superconducting side). 



The 'two fluid' picture here requires Eqs. (||) and (|^) 
be equal be satisfied in order to guarentee global current 
conservation. Once the junction geometry and applied 
voltage is specified, the only free parameter in Eqs. (||) 
and (M) is the superfluid velocity Vg — hq/m. Since 
the quasi-particlc transmission and reflection coefficients 
themselves depend on Vg, equating Eqs. (||) and (|^) is 
then a globally self-consistent procedure for determining 
the superfluid flow velocity Vg. Using this value for Vg, 
one then uses either Eq. (|^) or (||) to find the terminal 
currents for each value of the bias voltage V. 



B. Two-Fluid Approximation for Ballistic NS 
Junction 



We now restrict our attention to ballistic NS junc- 
tions, and approximate the energy bands near the Fermi 
level as simply rigidly shifted in energy by an amount 
±\A\{vs/vd), as determined from Eq. (|^). When the su- 
perconductor carries a finite supercurrcnt, we can then 
obtain the Andreev reflection coefficient from 



TnhME) = Ra{E - \A\{vg/vd)). 



(8) 



Yieie£^{E) is the Andreev reflection coefficient found by 
BTKllil when the superfluid flow is zero {vg = 0) , namely 



i?° (S) = <^ fE-^W-^ 



^E+^E^-\A\ 



\E\ < |A| 
\E\ > |A| 



(9) 



This rigid shift in the Andreev reflection coefficient, cor- 
responding to the rigid shift in the energy bands near the 
Fermi level, is shown in Fig. ^ Simply shifting the re- 
flection coefficients in energy is not a valid approxinmtion 
when a tunnel barrier is present at the NS intcrfaccEj, as 
also noted in Fig. ^. Since the differential conductance 
of an NS junction is 



dl 
dV 



^'^[l + RA{E^eV)~RN{E^eV)], (10) 
h 



a differential conductance measurement producing a 
larger than expected energy gap could point to signifi- 
cant superfluid flow in the junction. 

Shifting the Andreev reflection coefficient in energy 
makes it possible for the 'negative energy' Andreev re- 
flections, i.e. the Andreev reflection probabilities having 
_E < in Eqs. (^) to contribute to the excess electrical 
current. These 'negative energy' Andreev reflections are 
shown as the additional area under the Andreev reflec- 
tion probability for E > when > in Fig. |3| At zero 
temperature, Eqs. ([l]) and (||) give 

I = i2e/h)MN / [1 + K{E ~ \A\ivg /vd))]dE. (11) 
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To determine the superfluid velocity Vs 
Eqs. (^) and (0) require 



in Eq. (ITT]) 



\A\a{vs/vd) 



R\{E-\A\{vJvd))dE. (12) 



Equation ( |12[) is a self-consistent equation for the su- 
perfluid velocity Vs, and depends on the ratio of the 
number of conducting modes in the superconductor to 
the normal conductor a = Ms /Mm- The largest possi- 
ble excess current would occur if we could fix — s- oo 
in Eqs. (11), and would give twice the BTK result of 
7™°^ = 27*^^ = (8/3)(2eA//i). However, this theo- 
retical maximum excess current is not possible due to 
the constraint that the superfluid velocity must remain 
smaller than the depairing velocity. One must therefore 
discard any solution of Equation (n3) giving Vg > Vd- 




Energy (E/Aq) 

FIG. 3. Andreev reflection probability for an electron inci- 
dent on an NS interface when the superfluid velocity is zero 
(solid line) and when vs ~ Vd/2 (dashed line). Superfluid 
flow simply shifts the Andreev reflection probability by an 
amount /S.{vs/vd)- If a tunnel barrier is placed at the NS in- 
terface (dotted line), there is no simple relation between the 
Andreev reflection probabilities with and without superfluid 
flow. 



interface is formed where the narrow conductor meets 
the wide superconducting reservoir. Geometrical dilu- 
tion of the supercurrent where the superconductor widens 
into a thermodynamic reservoir moves the NS interface 
so that a stable point contact junction is formed. Forc- 
ing the narrower superconducting region into the normal 
state therefore creates a point contact NS junction hav- 
ing a = oo, i.e. the BTK limit of an NS junction. Forcing 
the narrower superconducting region to become normal 
therefore forces the excess current to fall abruptly to the 
BTK hmit shown in Fig. |(a). 
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C. Numerical Evaluation of Two-Fluid Formula 



We plot the numerical solution of Eqs. (|TT])-(|T^ for 
the total current through the NS junction as a function 
of the voltage in Fig. ^. When the transverse mode ra- 
tio is small, namely when a < 7/3 shown in Fig. ^(a), 
there is a voltage above which current conservation re- 
quires that Vs exceed the Landau depairing velocity of 
the superconductor. When Vg > Vd, the narrower super- 
conducting region (between the large normal contact and 
the large superconducting reservoir) becomes a normal 
conductor. This collapse of the order parameter in the 
narrower superconducting wire continues until a new NS 



FIG. 4. (a) When the transverse mode ratio a < (7/3), 
the order parameter in the narrow superconducting wire can 
collapse, giving rise to discontinuities in the I-V relation, (b) 
When Q > (7/3) no such discontinuities arise. The excess 
current is larger than the BTK result for a = 7/3, shown in 
both (a) and (b). The I-V evolves smoothly into the BTK 
result for a point contact NS junction when a oo. 

An excess current larger than the BTK limit is also 
shown in Fig. §(a), due to the additional 'negative en- 
ergy' Andreev reflections. When the transverse mode 
ratio is larger, namely when a > 7/3 shown in Fig. ^(b), 
the superfluid velocity in the narrower superconductor 
is always less than the Landau depairing velocity. Con- 
sequently, no abrupt drops in the current occur for any 
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value of the voltage. The excess current simply decreases 
gradually from its maximum value at a = 7/3 to the 
BTK value for a NS point contact at a = oo. 

We plot the maximum excess current in a ballistic NS 
junction versus the mode ratio a in Figure I The excess 
current at any given voltage is defined as the difference 
between the current carried by the NS superconducting 
junction and a normal NN junction, namely Iex{V) — 
lNs{y) ~ Inn{V)- The maximum excess current occurs 
at voltage for which I^x = Max [/ea;(F)]. For a < 4/3 
the excess current is given by the point contact value, 
namely the BTK result of (4/3) (2e A/ft,), and occurs at a 
voltage V — oo. When 4/3 < a < 7/3 the excess current 
occurs at a finite voltage < oo, and is larger than 
the BTK result. As the transverse mode ratio increased 
above a > 7/3, geometrical dilution of the supercurrent 
reduce the band tilting in the superconductor, reducing 
the maximum excess current of the junction. When the 
mode ratio is very large, so that a —^ oo we recover the 
point contact result for the excess current. 
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FIG. 5. The excess current Jcxc of a ballistic nor- 
mal-superconducting wire at zero degrees Kelvin is maximum 
when the transverse mode ratio is a = (Ms/Mn) = 7/3. 
/cxc approaches the point contact (BTK) value of 4/3(2eA/^) 
when Q — > 00. Collapse of the order parameter in the narrow 
superconducting wire limits 7cxc when a < 7/3. 



The circled dots in Figure 0|Show 



the results of a more 
realistic numerical|-Galculatior£3 for the excess current. 
As detailed in Ref.lla, we permit interband scattering at 
the NS junction and allow different conducting modes 
in the superconductor to carry different amounts of su- 
percurrentQ. The general behavior for the excess current 
versus mode ratio a of this more realistic NS junction 
model is quite similar to our simplified model, except the 
maximum in the excess current is slightly smaller and 
shifted to a slightly larger transverse mode ratio. The 
slightly lower excess current arises because the higher ly- 
ing modes have a smaller Fermi velocity, and therefore 
the Andreev reflection coefficients for the higher lying 
modes do not shift in energy as much as the lower ly- 



ing modes. We conclude that this more realistic model, 
even though still a two-fluid type model which only glob- 
ally conserves electrical current, confirms the essential 
features of our simpler independent mode calculation in 
this section. 



D. Limiting Cases of Two-Fluid Formula 



We can understand how the excess current lexc de- 
pends on the mode ratio a — Ms/AIjy in Figure |^ by 
examining Eqs. (|ll|)-(|l^) in different limits. The integral 
of the shifted Andreev reflection probability in Eqs. (O)- 
( |l^ ) can be done analytically to yield 



R''jE-\A\{vJvd))dE = eV 



(13) 



when eV < |A|(I + {vs/vd))- In this limit, Eqs. (|T|) and 
( p^ ) show the superfluid velocity increases linearly with 
bias voltage V. For larger biases, namely when eV > 
|A|(1 -|- (vs/vd)), the superfluid velocity increases more 
slowly with voltage, as determined from 



eV 



R\{E-\A\{vJvd))dE 



|A|(I + (t,./t;rf)) + |A| 



The factor 7 in Eq. mm is 



7 



cosh ( — 



2 ^6 ^ 



(14) 



(15) 



The excess current we obtain from 

hxc = (2e//i)Mjv / R\{E - \l\\{vslvd))dE. (16) 

Consider first the case where the narrower supercon- 
ducting wire is not driven normal, so that Vs < Vd for 
all values of voltage. In that case, the maximum excess 
current occurs when = 00, so that 7 = 00 in Eq. (p^). 
Eqs. ( p^ ) and (Q) for the superfluid velocity then reduce 
to (a — l)(vs/vd) = 4/3. The maximum allowed super- 
fluid velocity, (vg/vd) = I, then occurs for a transverse 
mode ratio of a = 7/3. The excess current from Eqs. ([l^ ) 
and dl3) then becomes 



2e|A| 



N 



a 



a ~ 1 



(17) 



when a > 7/3. The excess current reaches its maxi- 
mum value of lexc = (7/3)(2e|A|//i)MAr when a > 7/3 
as shown in Figure ^. Taking the limit a — > 00 in Eq. ( p^ ) 
recovers the BTK resuh lexc = (4/3)(2e| A|//i)MAr. 

We can also obtain an analytical solution for the ex- 
cess current using Eqs. (|l^) and (O) to determine the 
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superfluid velocity as {vs/vdj = eV/\A\a. The ex- 
cess current then follows from Eqs. (jl^) and ( |l3| ) as 
lexc = {2e\A\/h)MN\A\a{vs/vd)- Using the maximum 
alUowed depairing velocity of {vg/vd) — 1 then gives 

Ie.c = (^Af^) «■ (18) 

The range of allowed mode ratios a for which Eq. (|8|) is 
valid lie along the curve (vg/vd) = 1 — eV/\A\a. Further- 
more, Eq. ( p^ is only valid for ey/|A| < (1 + (vs/vd))- 
Combining all these requirements restricts the mode ra- 
tio between 1 < a < 2. However, for 1 < a < 4/3 the 
order parameter in the narrower superconducting wire 
collapes before the excess current reaches the BTK value. 
Eq. (^8|) therefore describes the excess current between 
4/3<Q;<2as shown in Figure |[ 

III. EXACT SOLUTION OF A SINGLE MODE NS 
JUNCTION 



In this section we wish to evaluate several assump- 
tions made in the two fluid model of section |l[ To 
do this, we solve the BdG equation, together with the 
self-consistency requirement for the order parameter, in 
a single band (ID) model of the NS junction where 
Mjv = Ms = 1 (a = 1). This self-consistent solution 
allows us to demonstrate how the supercurrent flow de- 
velops naturally from a self-consistent solution of the NS 
junction under a voltage bias. We use this solution to 
verify the approximate (two-fluid) procedure we use to 
guarentee global current conservation in section H One 
might expect that conserving the current only globally 
certainly gives qualitatively correct answers for the I- 
V characteristic. However, since the Andreev reflection 
probability does not vary much (as a function of energy) 
if we allow the order parameter to reach its final self- 
consistent form, the two-fluid model also gives accurate 
quantitative estimates for the I-V relation. 

Our self-consistent solution of the BdG equations in 
this section verifies the main assumptions used in our 
two- fluid model (when a = 1). When a voltage is ap- 
plied to the ballistic NS junction junction, the magni- 
tude of A(a;) remains approximately constant inside the 
superconductor (at zero temperature). However, the or- 
der parameter phase varies approximately linearly inside 
the superconducting metal. There are essentially no ad- 
ditional quasi-particles injected into the single moded NS 
junction at zero temperature, so the total current is sim- 
ply / = envs- The slope of the phase is related to the 
supercurrent velocity as d(j)/dx ^ 2q — 2vsm/h. The 
superfluid velocity is linearly related to the voltage as 
/ = {4:e^/h)V = envg. As voltage bias increases, the 
slope of the phase d(j)/dx increases until the superfluid 
velocity reaches the Landau depairing velocity, Vg = Vd 



(or d(j)/dx = l/'Co)- At this voltage V = Vc the order- 
ing parameter inside the narros superconductor collapses, 
and a stable point contact NS junction forms inside the 
wide superconducting reservoir. 

A. Self-Consistent Solution Procedure for BdG 
Equation 

The motion of quasi-particles in our one band NS junc- 
tion, including the superfluid flow inside in the super- 
conductor, is determiaed from the ID time independent 
Bogliobov-de Gennest^ (BdG) equation 

{H{x)-^i A{x) \ f u{x)\ _ f u{x)\ 

[ A*{x) -{H*{x)-^,}J[vix)J-''[v{x)J 

(19) 

The one-electron Hamiltonian H{x) in Eq. ([l^ ) is 

Hix)^~^-f^ + V{x). (20) 
2m dx'' 

The order parameter A(x) in Eq. ( [T9| ) is given by 
A(a;) = 9{x)F{x) 

= 9{x)^Vpn{x)Upn{x)fpnO{\Epn\ -TiUJd), (21) 
pn 

where g{x) is the electron-phonon interaction strength at 
each point and F{x) is the pair correlation function. (Al- 
though we assume g{x) is local, in reality it is spread over 
a correlation distance vp/i^D-) The index p in Eq. ( [2l| ) 
denotes the lead from which the scattering state origi- 
nates, namely p = TV is the left lead and p = S the right 
lead. The index n in Eq. (^) denotes the good quantum 
numbers in the lead, namely n = (fc, (3) where k is the 
wavenumber and (3 = (e, h) the electron-like and hole- 
like states. The sum in Eq. ( ^l]) runs over states injected 
from the leads, including both positive (En > 0) and 
negative energies {En < 0). The coherence factors u{x) 
and v{x) in Eq. (UW are functions of the order parameter 
A(x) through Eq. 

In this section we show the self-consistent solutions of 
Eq. (|l^) and ( |2l| ) for a voltage-biased NS junction. De- 
tailfipf the self-consistent solution procedure are given in 
Ref.a To solve the order parameter self-consistently, we 
flrst assure an initial or zeroth order guess Ao(a;) for the 
order parameter. We then divide the one dimensional 
space is into differential elements, where the magnitude 
of the order parameter superfluid velocity are constant 
in each section. We match the wavefunctions and their 
derivatives at each interface to obtain the zeroth order 
wavefunctions uq{x) and vo{x). The flrst iteration for the 
order parameter Ai(x) we then obtain from Eq. IGV) us- 
ing the zeroth order wavefunctions uo{x) and vq{x), etc. 
The zeroth order guess for the order parameter Ao(a;) 
can either be constant, i.e. Ao(a;) = A, or it can contain 
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a superfluid flow, i.e. Ao(a;) = Ae^"'^. Given the same 
electrical current flow, either initial guess for the order 
parameter converges to the same final answer. 

The voltage bias across the NS junction is a boundary 
condition which determines the Fermi occupation proba- 
bilities fpn in Eq. ( pT| ) . The occupation probability /p„ of 
a scattering state (p, n) which originates inside the nor- 
mal or superconducting reservoir p, is the same for holes 
and electrons when the applied bias is zero {V = 0). 
Under a voltage bias, however, electrons in the normal 
metal are occupied up to an energy /ip+i eV , while holes 
are occupied up to an energy ^ — eVs^ These different 
Fermi factors electron-like and hole-like quasi-particles 
injected from the normal metal are shown schematically 
in Fig. H The unequal occupation probabilities for holes 
and electrons injected from the normal metal causes these 
two classes of scattering states to contribute differently 
to the sum in Eq. ( ^ ) under an applied bias. We can 
write the Fermi factors as 



fpn = fpp = f{E - eVpp) , 



(22) 



where ,f{E) = + exp (£;//csr)]. Here eVpp is effec- 
tive biasing voltage (or effective electrochemical poten- 
tial) applied to the (p/3)th lead, namelytj 



VNe = V 



VNh 



-V . 



(23) 
(24) 



In this paper the superconducting leads are grounded so 
that Vsp = 0. 

To obtain the electrical current /(x), we do not in- 
voke any ad hoc 'source term' as done in RefO, but in- 
stead simply evaluate the electrical current operator for 
a scattering state originating in lead q (having quantum 
number n) and terminating in lead p, namely 



E 

pn 



{Ju + Jv)a:r)n fp 



E 

pn 



(25) 



The Ju and are Schrodinger currents associated with 
the waves u and v, namely J„ = (e?i/m)Im{u*(a;)VM(a;)} 
and J.„ = (eh/m)l'ni{v*(x)yv{x)}. The 'vacwum current' 
due to the filled hole band is argued in Ref.llHl to be zero, 
namely X]gn("^y)p:9n ~ have confirmed for the 



NS junction. Solving Eq. (gj) togethcrpWiiJ|--Eq. (|21 
guarentees electrical current conservationl3ij,Ll,llIllij, even 
when the superconductor is far from equilibrium. A proof 
of this statement for NS junctions follows from generaliz- 
ing the discussion in Appendix B of RefB to the nonequi- 
librium case. If the nonequilibriuni-aystem involves two 
superconductors at different biasestS, current conserva- 
tion is more complex and J2qni'^'u)p;qn 0. 



B. Order Parameter Phase 

Consider first the NS junction, where the coupling con- 
stant ^(a;) is 



9r x>0 
x<0 



(26) 



We choose and so that the critical temperature of 
the right superconductor is Tc — 6.6K. Our initial guess 
for the order parameter we take to be 



Ao(x) 



A x>Q 
X <0 



(27) 



We therefore do not force a superfiuid flow inside the 
superconductor from our zeroth order guess for the order 
parameter. 
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FIG. 6. Both (a) the electrical current throughout an NS 
junction, and (b) the phase in the superconducting metal, 
develop naturally from a self-consistent model. The ini- 
tial order parameter guess Ao(a;) assumed zero supercurrent, 
however a superfluid flow appeared naturally upon reaching 
self-consistency. 

Fig. ^(a) shows the electrical current l(x) versus po- 
sition X inside the superconductor. The numbers beside 
the lines in Fig. ^(a) denote the iteration number. For 
the first iteration [N = 1), the electrical current dies 
off within a coherence length of the NS interface, so that 
the electrical current is not conserved. After the iterative 
scheme converges {N = 700) in Fig.|(a), we see the elec- 
trical current I{x) is constant as a function of position. 
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indicating the electrical current is indeed conserved. In 
Fig. 11(b) we plot the order parameter phase inside the su- 
perconductor as a function of position. A uniform phase 
gradient develops inside the superconductor when the it- 
erative scheme has converged to self-consistency, showing 
that the development of a supercurrent is necessary to 
guarentee current conservation. We expected this con- 
stant order parameter phase gradent in the NS junction, 
since it is similar to that constant order parameter phasa 
found in the self-consistent solution of the SNS junction.EI 
We can understand these result using our two-fluid pic- 
ture from section and assuming rigidly shifted energy 
bands. The maximum current of 4eA//i is reached when 
the voltage eV — A. At this point the energy bands 
have shifted up faster than the Fermi level of the nor- 
mal contact, and thus no direct transmission of electrons 
across the junction inside the energy range < -B < eV 
is allowed. Andreev transmission will also be small, as is 
usually the case in NS junctions. There will therefore be 
essentially no quasi-particles above the Fermi level of the 
superconductor, while all states below the Fermi level 
are filled. The quasi-particle contribution to the cur- 
rent inside the superconductor, Iqp in Eq. (|^), is essen- 
tially zero in the single mode NS junction when eV < A. 
The electrical current is therefore / = envg in this single 
moded NS junction, the same as for a uniform ID su- 
perconductor. The order parameter magnitude collapses 
when the superfluid velocity equals the dcpairing velocity 
Vs = I'd, at a bias voltage eV = A. 

C. Order Parameter Magnitude 



Fig. ^a) shows the magnitude of the condensation am- 
plitude F{x) as a function of position in both the normal 
and superconducting metal at zero temperature. The 
solid line indicates a bias voltage oi V = 0, while the 
dashed fine is for a bias voltage V = 0.95T4 (^s = 0.95ud). 
The general form of the condensation amplitude F{x) for 
a ballistic NS junction at equilibrmrnjs well known from 
earlier non-self-consistent models.E2l,til. We find substan- 
tial agreement between these earlier results and our fully 
self-consistent calculations. In the superconductor, F{x) 
is suppressed from its bulk value near the NS interface. 
In the normal metal, F(x) shows behavior quite similar 
to the low temperature experimental results of MotaE3. 
The dotted line in Fig. ^a) ia i^(a;) oc l/{\x\ + xq), 
found experimentally by Motac3. The value of xo used 
in Fig. 0(a) is xq = ^o- The fit between the experimen- 
tal determined form F{x) oc l/(|x| -I- xq) and the results 
of our self-consistent calculation is quite good. The re- 
sult F{x) oc 1/(1 a; I + a^o) for the pair correlation func- 
tion at low temperature in an NS junction at equilibrium 
{V = 0) was also pointed out by FalkEil for the asymp- 
totic limit X oo. 




FIG. 7. The magnitude of the coherence function F{x) (a) 
changes little when a large flow is present in a NS junction 
with an applied voltage. The solid line is F{x) when the ap- 
plied voltage is zero. When the applied voltage is large .95A 
F{x) shows little change (dashed line) from the junction in 
equilibrium. [The phase of F{x) changes linearly with pos- 
tion throughout the NS junction.] (b) In an S'S junction at 
a temperature above the critical temperature of S', the order 
parameter is non-zero in the "normal" metal. The finite tem- 
perature in combination with a moderate supercurrent causes 
a supression of the larger gap superconductor, (dashed line) 

In an NS junction, the ordering parameter A(a;) van- 
ishes in the normal metal because the electron-phonon 
coupling constant g(x) is zero there. For a finite order- 
ing parameter A(a::) to exist inside the normal metal we 
must have g{x) 7^ in the normal metal. One way to 
achieve a non-zero g{x) in the normal metal is to fabri- 
cate an S'S junction, where the superconductor S' has 
a smaller critical temperature than S. If we then ele- 
vate the temperature so that Tc > T > T^,, we effec- 
tively form an NS junction where g{x) is not zero in- 
side the normal metal. Fig. ^(b) shows a self-consistent 
calculation for an such S'S junction, where Tc = 6.6K, 
= OMK, and T = 2K. Unfike the NS junction, where 
A(a:;) = inside the normal metal, we see a non-zero 
'tail' of the ordering parameter extending into the normal 
metal in Fig. 0(b). At zero temperature, the bulk value 
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of the order parameter inside the weaker superconductor 
is A'(r = 0) = 0.1A(T = 0). From Fig. @(b) we see that 
A(a; = 0, T = 2K) ~ 2.5A'(T = 0), larger than even the 
bulk value of the order parameter in the weaker supercon- 
ductor at zero temperature. When a voltage is applied 
to the S' S junction, namely eV = .7Ao in Fig. 0(b), the 
ordering parameter inside S is now suppressed from its 
bulk value at T = 2K. This degradation of the order 
parameter at finite temperature, when the superconduc- 
tor carries a finite supercturrent, is similar to that of a 
bulk superconducting wireEI. The tail of the ordering pa- 
rameter extending into S' is only slightly changed in the 
presence of the supercurrent. 



D. Local Density of States 



In addition to the magnetic susceptibility techniques 
used by Mota, which explore the condensation amplitude 
F{x) in the normal metal, another method to experi- 
mentally investigating how the ordering parameter A(a;) 
varies near NS mterfaces is tunnelling spectroscopy us- 
ing an STM tip.E^. We expect that a measurement of the 
differential conductance dl/dV at the STM tip is propor- 
tional to to the local density of states N{x,E). We can 
calculate the local density of states using the equation 



N{x,E) ^-Y\up,n{x,E)\ 



dk 



(28) 



where p is the lead index and the quantum number 
n = (fc, /?). Fig. |^(a) shows the local density of states for 
an NS junction having an applied voltage of eV — .TAq. 
At this bias voltage, the superconductor carries a super- 
current of Vs — 0-7vd- The solid line shows N{x, E) at a 
position X = 5^o inside the superconductor. The original 
peak in the density of states, which occurs at = A 
when the superfluid flow Vg = 0, splits into two separate 
peaks. As the energy bands inside the superconductor 
tilt under the superfluid flow, the band edges move to 
the energies E = [\ ± (ti<;/fd)]A, as do the peaks in the 
density of states. In the normal metal, the density of 
states is approximately constant for energies of interest. 
The constant density of states in the normal metal is due 
to a zero pairing potential A(a::) = inside the normal 
metal, since g{x) = in the normal metal. 



03 

Z) 

B 
S 

< 

d 
d 



- 










Superconductor 












^Normal Metal 




1 1 1 1 




c 

Z) 

>. 
^_ 

S 
< 

C/j 

d 



eV/Ao 

FIG. 8. (a) When a supercurrent is present in an NS junc- 
tion, the local density of states inside the superconductor 
shows two peaks. The local density of states inside the normal 
conductor is constant for all a; < 0, indicating the lack of any 
energy gap due to the proximity effect, (b) In an S'S junction 
when a supercurrent is present, where S' is a superconductor 
above its transition temperature, structure is seen in the local 
density of states inside the "normal" S' metal. 

Fig. ||(b) shows the local density of states for an S' S 
junction having Tc > T > T^. We evaluate the local den- 
sity of states at a position x — — in ttiG normal metal. 
The local density of states inside the stronger supercon- 
ductor is approximately the same as Fig. ^(a). For the 
two different applied voltages, where Vg = (solid) and 
Vs = .5vd (dotted), the presence of the superconductor 
changes the density of states inside the normal metal. For 
Vs ^ 0, there is a depression in the density of states near 
E ^ 0, showing the partial development of an energy 
gap inside the normal metal. The density of states does 
not go to zero at _B = 0, since quasi-particles incident 
from the left contact can still propagate to the position 
X = As the current increases, the density of states 
inside the normal metal becomes flatter due to injection 
of quasi-particles from the tilted energy bands inside the 
superconductor. The two small peaks at eV = .5Ao and 
eV = 1.5Ao when Vg = -Bvd are again associated with the 
tilted energy bands inside S, and are significantly broad- 
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ened by thermal smearing at T = 2K (which we have 
ignored in Fig. ||(b)). 

To summarize, while looking to justify our two-fluid 
model for the effect of superfluid flow in an NS junction, 
we studied the pair corrlation function F{x) and order- 
ing parameter A{x). We looked at F{x) and A(x) inside 
both (1) an NS junction at temperature T = 0, and (2) 
an S'S junction at a temperature Tc > T > T^. Here S' 
is a superconductor having an order parameter smaller 
than S, so that < Tc. The S'S junction at this tem- 
perature is therefore a type of NS junction. First, the 
condensation amplitude is F{x) ~ xo/{\x\ +xo), approx- 
imately independent of the applied voltage. Second, at 
the voltage Vc ~ Ao/e, the ordering parameter of the 
superconductor collapses, i.e. Aq — > 0. The supercurrent 
carried inside the superconductor at a voltage V = Vc 
is approximately equal to the Landau depairing current 
Ic = envd- Third, we show how the supercurrent changes 
the 1-D local density of states per unit energy N{x, E) at 
various points in the normal and superconducting metals, 
both for the NS junction and the S'S junction. The lo- 
cal density of states shows the influence of the superfluid 
flow. 



E. Locally (But Not Globally) Gapless 
Superconductivity 

The phase gradient of F{x) is approximately constant 
for both the NS junction in Fig. |^(a) and the S' S junction 
in Fig. R(b). Since the energy gap A (a;) decays to zero 
in Fig. CTa)-(b), there are regions of local 'gapless' super- 
conductivity near the NS interface where Vs > A{x)/pF. 
However, these regions of local 'gapless' superconduc- 
tivity do not affect bulk properties such as critical cur- 
rent, critical temperature, etc. Since our self-consistent 
model shows that the supercurrent is approximately con- 
stant throughout the superconductor, depairing will oc- 
cur when Vs = Vd everywhere in the superconducting 
wire. 

A different type of 'global' gapless supfirconductivity 
has been discussed by Sanchez and Solso. In this pro- 
posed gapless superconductorQ, a superconducting wire 
in contact with an NS junction is postulated to exist 
when Vg > Vd- It is true that such a non-equilibrium self- 
consistent solution to the BdG equations with v^ > Vd 
does exist for a uniform superconducting wire. However, 
connecting such a wire in this novel 'global' gapless su- 
perconducting state to an NS junction imposes the ad- 
ditional constraint that the self-consistency condition in 
Eq. (21) must be satisfied at every point in space. When 
Eq. (21) is satisfied, the magnitude of the order parame- 
ter can no longer be a constant in space, as required for 
the globally gapless solution proposed in Ref.cl to exist .j^ 

Another way to view the situation proposed in Ref.cl 
is that, for a given applied voltage, the constraint of 
current conservation fixes both A(x) and Vs, leaving no 



more degrees of freedom. The order parameter A(a;) 
and superfluid velocity Vg cannot be adjustedjndepen- 
dently as required for the bulk solution of Ref.Q to exist 
in an NS junction. Our self-consistent model shows in- 
stead that the ordering parameter of the superconducting 
wire collapses when the junction voltage is approximately 
eV = A, or equivalently when v^ = Vd- In short, we see 
no possible way to achieve the non-equlibrium conditions 
necessary for this novel gapless superconducting state of 
Ref.Q to exist by connecting a superconducting wire to 
a ballistic NS junction. Fortunately, many other inter- 
esting measurements are possible at ballistic NS inter- 
faces without invoking the global gapless superconduct- 
ing state of Reffl. In particular, the excess current larger 
than the point contact limit, which we find in section [ 
in no way depends on the gapless state proposed in Ref 



IV. CONCLUSIONS 



It is possible to experimentally observe excess currents 
larger than the (4/3)(2eA//i) found for ballistic NS point 
contacts. By varying the width Ws of a superconducting 
wire in contact with a normal metal having width Wn, 
one can vary the effect of the supercurrent on the energy 
bands in the superconductor. Varying the widths of the 
two conductors controls the ratio of the number of con- 
ducting modes a ~ Ws /Wn in the ballistic NS junction. 
We find the excess current attains a theoretical maxi- 
mum of (7/3)(2eA//i) when a ~ 7/3. For 1 < a < 7/3 
it should be possible to observe discontinuities in the I- 
V relation of the NS junction when the superfluid ve- 
locity Vs exceeds the Landau depairing velocity Vd- Al- 
though these results follow from a simple model which 
treats all conducting modes as equivalent, we confirmed 
the qualitative results using a more realistic model which 
includes the different supercurrent carried in each con- 
ducting mode and the scattering between the different 
modes. These predictions are based on a 'two-fiuid' solu- 
tion of the BdG equations, in which current conservation 
is violated locally near the NS interface. 

To confirm our that our 'two-fluid' type treatment of 
the superfluid flow in NS junctions generates qualita- 
tively accurate predictions for the I-V relations, we solved 
the BdG equations self-consistently for a single mode 
NS junction under an applied bias. Current conserva- 
tion follows automatically in this self-consistent scheme, 
and shows that the superfluid velocity is indeed constant 
throughout the NS junction. The two important features 
confirmed in this self-consistent solution are (1) the su- 
perfluid velocity and terminal currents are the same as 
required by our 'two-fluid' scheme and (2) the order pa- 
rameter indeed collapses when the superfluid velocity 
equals the depairing velocity Vd- We did not perform 
a completely self-consistent calculation for the multiple 
moded NS junction, as we believe all the essential ele- 
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merits of this problem (so far as the excess current is 
concerned) are encompassed in the two-fluid treatment 
of the multiple-moded NS junction. 

Having obtained the a self-consistent solution of the 
BdG equations for an NS junction under bias also enabled 
us to study the pair correlation fmiction F{x), order pa- 
rameter A{x), and local density of states N{x, E) in the 
NS junction. The peak near the superconducting energy 
gap at -E = A in local density of states iV(a;, E) inside 
the superconductor is split by the superfluid flow, as can 
be measured using STM spectroscopy. If the electron- 
phonon interaction in the normal metal is nearly zero, 
the density of states in the normal metal is unaffected 
by the presence of the superconductor. If the normal 
metal iV is a weak superconductor held above its tran- 
sition temperature, the STM can also measure changes 
in the local density of states N{x, E) inside the normal 
metal. Due to the uniform superfluid velocity inside an 
NS junction, and the reduction of the order parameter 
near the NS interface, the Landau depairing condition is 
locally violated and a type of gapless superconductivity 
occurs locally near the NS interface. 



The derivation of the electrical current in RefJi ap- 
plies for an NS junction subject to a superfluid flow. 
EvalxLating the current operator inside the normal region, 
Rcf.tj finds the well known formula developed by BTKO, 
namely 

/oo 
[l-TNe,Ne{E)+TNhMe{E)] 
-OO 

[f{E - eV) - f{E)]dE. (29) 

Here we use the notation of Rcf.Ei, where Tn^.nc — Rn 
and Tjyh.Ne — Ra are the particle current reflection prob- 
abilities for an electron-like quasi-particle incident from 
the normal metal (right index, We') to transmit as a 
hole or electron in the normal metal (left index). The 
sum rule from Eq. (^ then gives / = Iqp + I a, with 



IqP = {2e/h) I [Tse.NeiE) ~ Tsh.NeiE)] 

[f{E - eV) - f{E)]dE, 



(30) 



and 
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To evaluate the electrical current operator inside the 
superconductor, we first note that the energy bands, in- 
side superconductor subject to a superfluid flow area 



2m 



2„2 



2m 



/i = ±J{ETfi' kq/my 



(32) 



VI. APPENDIX: DETERMINING THE 
SUPERFLUID FLOW VELOCITY 

If the pairing potential and electron wavefunctions are 
determined from a self-consitent solution of the BdG 
equations, (1) electrical current will be conserved every- 
where in space and (2) the superfluid flow will develop 
naturally as the scheme evolves towards self-consistency 
(c.f. section 



III). In this appendix we examine a differ- 



ent 'globally self-consistent' scheme, where current con- 
servation is only guarenteed at the device leads. In 
this scheme, the correct value of the superfluid veloc- 
ity vs = fiq/m is determined by equating the current 
operator evaluated on the normal side of the NS junction 
with the same operator evaluated deep (several coher- 
ence lengths) inside the superconductor. We then adjust 
the superfluid velocity (which is a free parameter in the 
scheme) until the current flowing out of the supercon- 
ducting lead is the same as the current flow in the nor- 
mal lead. We consider only a single moded NS junction 
with Mjv = Ms = 1, although the results here are easily 
generalized to multiple conducting modes. 



The discussion in Ref.S can be extended to show that the 
particle current incident from the superconductor Jp — 
(\/Ti){dE / dk). Thus, the particle current incident from 
the superconductor per unit energy is simply \/h. Quasi- 
particles incident from the superconductor will then carry 
an electrical current of the form 



dE 



Tout Jon 

•Jp ■-'q 



T" 



rout 
Jp 



(33) 



We recognize the particle current transmission coefficient 
in Eq. (p3[) as 



out, in ^ 



(34) 



To evaluate the Jq^^LJp^^ term in Eq. (|3^) we use sev- 
eral results from RefJS. The scattering states inside the 
superconductor have solutions of the form 



u{x) 
v{x) 



f Ukq exp(zga;) \ ^^kx 
\vkq exp{~iqx) 



(35) 
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The electrical current carried by each occupied state is therefore 

M)f, (36) 



and the particle current for each occupied state is 

Jp = (^) {\uk,\' - \v,,n f + {\u,,\' + \v,,\') /. (37) 

The state are normalized so that \ukq\'^ + \vkq\'^ ~ 1- Hence the factor 

1 ^ (g/fc)(K,p-Kn 



(38) 



The first term in Eq. (^8|) one can show is simply the ratio of the density of states in the superconductor to that of 
the normal metal (for a fixed value of wavevector k, not a fixed energy), namely 

^'^^^ = ^ (|Ufc,|2-|t,,,|2) + (g/fc)- (^^^ 

We make the translation between wavevector k and energy E inside the superconductor using Eq. (|3^). The second 
term in Eq. ( |3^ ) will be only a minor correction to the first, being nearly zero near the Fermi level, equal to (g/fc) << 1 
over most of the energy range, and equal to 1 only near the bottom of the electron energy bands. The first term in 
Eq. (^) therefore dominates, being much larger than 1 near the Fermi level and equal to 1 over most of the energy 
range. I_. 

Applying the procedure outlined in Ref.li3 for construction of the scattering states, and multiplying by their appro- 
priate Fermi occupation factors, we find a total current inside the superconductor of 

Is^h+h- (40) 
The current Ii arises from first term in Eq. ( |3^ ) and is 
2e 

^1 = 7^ y_ [Nl^,^,iE)TseME) - Nl^^,^^{E)TsHMmf{E - eV) - f{E)]dE 

Op 2e f°° 

+ NLoutmV{E)^l\dE-- Nl,^„,{E)[2f{E)^l]dE. (41) 

To obtain Eq. (41) we used the sum rule and electron hole symmetry, Eqs. (C.l) and (C.6) of Ref.El. We distinguish 
between Ng,^ ^^^{E) and Ng^^^{E) in Eq. (p|), since the incoming and outgoing electron- like quasi-particles have 
different densities of states. The second term in Eq. ( p8| ) results in a small correction current I2, proportional to 
{q/kp). The rather cumbersome expression for I2 is 



_2e 

2e 
~ T 
2e 

_ ^ 

We again translate between k and E inside the superconductor using Eq. ( |32| ) , taking care to assign the appropriate 
branch of the dispersion curve for incoming or outgoing electron- or hole- like particles. The first two terms in I2 are 
a small correction to the superfluid flow due to additional quasi-particle injection, while the last two terms are a small 
correction to the equilibrium superfluid flow. A rigorous treatment guarenteeing global current conservation at any 
temperature would equate the current / from Eq. ( p9| ) with the current Is from Eq. (^^, adjusting the superfluid 
velocity q until I = Is- 



[{q/k){\Ukg\^ - \vkqn]se,outNle,outiE)Tse,NeiE)[fiE - eV) ~ f{E)]dE 
[{q/k){\ukg\' - \vkqf)]sh,outNlf,,outiE)Tsh.Ne{E)][f{E - eV) - f{E)]dE 
[(g/fc)(|ufc,|2 - \vuq\')]se,outNl,out{Em{E) - l]dE 

[(g/fc)(kfc,P - \Vkq?)]se,rnNl,,AE)[2f{E) ~ l]dE . (42) 
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We now analyze the validity of our 'two-fluid" proce- 
dure from Section We henceforth neglect the current 
I2 as insignificant compared with Ii. At zero tempera- 
ture, the current operator evaluated inside the supercon- 
ductor we find from from Eq. (^) as 



2e 



eV 



/l (T = 0) = ^ / [7V|,,„„, {E)Tse,Ne [E) 



Nl^,out{E)Tsh,Ne{E)]dE 
4eA 



(43) 



The zero temperature limit of Eq. (^9|) is given in Sec- 
tion H as Eqs. The second term in Eqs.(|3|) is 
the superfluid flow term Ic- So we can certainly identify 
I A from Eq. with Ic- However, Eq. (0) for Iqp is 
not exactly equal to the first term in Eq. (Esh, the differ- 
ence being the additional factor of the superconducting 
density of states N^{E) in Eq. Q. 

For the Cooper pair flow away from the NS interface 
(which we are considering in this paper), and for the 
energy range between and eV, the outgoing hole-like 
quasi-particle conduction channel opens at a lower energy 
than the electron- like quasi-particle channel (see Fig. g). 
This means the flrst term in Eq. (BS) will be larger and 
more negative than Iqp from Eq. IW, requiring a larger 
value of the superfluid velocity Vg at each value of the 

f plied voltage V than in the two-fluid model of Section 
We conclude that the treatment in Section || therefore 
underestimates the effect of superfluid flow on the excess 
current. 
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